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Q\ , Abstract 

p ■ It is noticed that the excellent proof of the connection of magnetic flux 

'^ ! quantization and off-diagonal long range order(ODLRO) presented recently 

o 
o 



X 



by Nieh, Su and Zhao suffers from an imperfection, namely, the f-factors in 
the case of finite translation do not satisfy /(a)/(b) = /(a + b), which was 
imployed in the proof. A corrected proof is proposed to remedy this point. 
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The study of ODLRO of fermions, especially those strongly correlated, has drawn much attentic 
recently [in -[|1^ since the onset of ODLRO implies a phase transition in the system. This trend 



investigations is currently believed to be able to sheed some light on the explanation of the high-! 



superconductivity in cuprates by a number of physicists. In his original paper |T3], Yang argued th, 



the existence of ODLRO will likely lead to Meisner effect and magnetic flux quantization. The k( 
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spirit of the argument asserts that the dependence of the thermodynamic partition function on tl 
flux, which provides phase factors responsible for space tranlations of the reduced density matric 
, is very similar to the dependence of Bloch wave function on the wave vector in crystals. Althouj 
this argument is intuitive, it is not clearcut and direct. Quite recently, Sewell [p^] gave a proof of tl 
relationship of Meisner effect and ODLRO. Later, Nieh, Su and Zhao |jT9| extended the proof to tl 
case of flux quantization. Yet, their extention suffers from an imperfection when the lattice potenti 
is taken into consideration. This note is to remedy this point. 



In general, the reduced density matrix P2(i'i, ^^2] ^i, r2) can enjoy a spectral seperation||T^ 



P2(r'i, r'2; ri, rs) = J2 "^"^^(ri, ^'2)^1(^1, ra) ( 

s 

where $s(ri,r2) are normalized eigenfunctions of p2 with eigenvalues Og'- 

J ^I{ri,r2)^s{ri,r2)dridr2 = 1 ( 

/ /92(ri, r'a; ri, r2)dridr2$s(ri, rg) = a^$s(ri, r^) ( 

The existence of ODLRO in p2 means that in the off-diagonal long range (ODLR) limit, i.e., 

r'i — Yj |— > 00 while keeping | ri — r2 | and | r[ — r'g | finite, p2 behaves as 

P2{r[, r'2; ri, r2) ^ a^{r[, r'2)<l>*(ri, r2) ( 



where a is the largest eigenvalue which is of order 0{N). The important observation made in[p!9| 



that the space translation r ^ r — a will effect a gauge transformation, based on which two comple 
sets of basis functions were obtained. As a consequence, the transform property of p2 is 

ic 
P2{r'i, r'^] Tl, ¥2) = exp{ — [xa(ri) + Xe^i^i) - Xa(ri) - Xa(r2)]}p2(ri - a, r'^ - a; ri - a, r2 - a) (. 

where Xalr) = a • Ao(r) + 0(r — a) — (t>{r), Aq = |B x r, 0(r) is a gauge function which is in gener 
multivalued. Therefore, in the ODLR limit, we have 

$(r;,r^)$*(ri,r2)^exp{ — [xa(r;) + Xa(r'2)-Xa(ri)-Xa(r2)]}$(r;-a,r^-a)$*(ri-a,r2-a) (1 



which implies 



t€ 

<l>(ri, r2) = /(a) exp{ — [xa(ri) + Xa(r2)]}$(ri - a, r2 - a) 



where /)a) is a position-independent phase factor, /(O) = 1. So $ differ by a phase factor at differe: 
space points. 

Let us consider any two finite successive space translations , firstly by a and then by b. Since 

Xb(r - a) = b • Ao(r - a) + (/)(r - a - b) - 0(r - a) (. 

we have 

Xa(r) + Xb(r - a) = (a + b) ■ Ao(r) + 0(r - a - b) - (p{r) - -b ■ B x a = Xa+b(r) - -b ■ B x a (' 
Accordingly, 

$(ri, rs) = /(a)/(b) exp{ — [xa+b(ri) + Xa+b(r2)]} exp{[— b ■ a x B]}$(ri - a- b, ra - a- b) (H 

Thus for a closed curve C, which is in general zigzag when the translation vectors a, b, ... are finit 
a + b + ... = 0, eq.(fO) together with the singlevaluedness of $, which is ensured by the singl 
valuedness of p2, give 

/(a)/(b)...exp{£B-S} = l (1 

where S is the area enclosed by C with the orientation of the translation vectors. It is obvious th; 

/(a)/(-a) = 1 (l: 

Now suppose that the three primitive vectors of the crystal are aj, 2 = 1, 2, 3. Consider now that tl 
closed path constituted by aj, a^-. — aj and — a^, we have then 

— 2?> 
/(a,)/(a,)/(-a,)/(-a,) = exp{^— B ■ S} (L 

The quantization of flux follows then from (12) 

B ■ S = n L 

2e ^ 



Since the size of such a region enclosed by any two primitive vetors is of the typical order (lA) 



2 



n ^ 0, B must be of the order of magnitude B ~ {lA}^'^^^ ~ 10^ G (Note that B is usually parall 
to S in experiments), which is much stronger than the typical critical field strength of less than li 
G, therefore n must be zero, hence B = 0. This states that if S is simply connected, i.e., the regie 
enclosed by C is superconducting everywhere, we have Meisner effect. On the other hand, if S 
multiconnected, i.e., it contains normal regions which are of mesoscopic sizes, n may be non-ze: 



integers. In this case, we have the quantization of magnetic flux. 

We can also draw the same conclusion in another way. Since (10) can be also expressed as 

$(ri, rs) = /(a + b) exp{ — [xa+b(ri) + Xa+b(r2)]}$(ri - a - b, rs - a - b) 



we have 



and similarly 



IP 
/(a + b) = exp{[-B • (b X a)]}/(a)/(b) (H 



IP 

/(b + a) = exp{[-B • (a x b)]}/(b)/(a) 



Since /(a + b) = /(b + a) and /(a) are c- numbers, we infer that 

,2ze^ , , , 
expj -;^B ■ (a X b) = 1 
Tic 

Therefore, for any two lattice vectors a and b, we always have 

■r-. / IN 2tiTic 
B ■ (a X b) = n- 



2e 



The imperfection of the proof in[|19| is that, when the lattice potentials are taken into accour 
the translation vectors can only be lattice vectors. But for lattice vectors a and b, the multiplicatic 
rule is (16) instead of the following 

/(a + b) = /(a)/(b) (21 

which was used in [|l9l as the basis to prove the flux quantization. Eq(20) holds only for infinitesim 
translations if neglecting second order infinitesimals. 

The relation is somewhat similar to that of magnetic translation operators which realize a r; 



representation of the translation group pO|. Consider the Hamiltonian of A^ electrons 



^ = Ei[P^ + >(r.)f + ^Enr.)+E^c(r.)+E^«^-B (2 

«=1 i,j i i 

Define j3 = ^ . For any lattice vector R^ = ■niiai + r?T,2a2 + rrasas, consider the operator 

^ i ^ r- 

r(R„) = exp{5] -R™ ■ P,} exp{z ^(/3 x R^) ■ ^} (2: 

We have TT^ = I and[T(R^),H] = 0. Since 

r(R„)T(R„) = exp{t^Np ■ (R„ x R„)}T(R„ + R„) (2, 



we have 

T(R„)T(R„) = exp{tNP ■ (R„ x R„)}T(R„)T(R„) (2^ 

This is what (16) is similar to. As the Bloch wave functions can be designated by the eigenvalues 
the translation operators, i.e., the wave vectors k and the Hamiltonian, in the presence of magnet 
field, the wave functions can be characterized by the eigenvalues of the Hamiltonian and the magnet 
translation operators belonging to any chosen commutable set TZ = {R | [T(Rm), T(R„)] = 0}. 
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